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in ■ ABSTRACT 

. We present results from a search for additional transiting planets in 24 systems 

already known to contain a transiting planet. We model the transits due to the known 
planet in each system and subtract these models from lightcurves obtained with the 
^ ■ SuperWASP survey instruments. These residual lightcurves are then searched for ev- 

' idence of additional periodic transit events. Although we do not find any evidence 

for additional planets in any of the planetary systems studied, we are able to charac- 
Tlj" ■ terise our ability to find such planets by means of Monte Carlo simulations. Artificially 

CO I generated transit signals corresponding to planets with a range of sizes and orbital 

\^ . periods were injected into the SuperWASP photometry and the resulting lightcurves 

■ searched for planets. As a result, the detection efficiency as a function of both the 
Q^ \ radius and orbital period of any second planet, is calculated. We determine that there 

■ is a good (> 50 per cent) chance of detecting additional, Saturn-sized planets in P ~ 

10 d orbits around planet-hosting stars that have several seasons of SuperWASP pho- 
tometry. Additionally, we confirm previous evidence of the rotational stellar variability 

' of WASP-10, and refine the period of rotation. We find that the period of the rota- 

5-H ' tion is 11.91 ± 0.05 d, and the false alarm probability for this period is extremely low 

^ ■ 10-13). 

Key vifords: 

planetary systems, techniques: photometric 



1 INTRODUCTION 

1.1 Multiple planet systems 

Of the 347 presently known extra-solar planets, 90 are 
known to reside within multiple planet systemfl All of these 
systems, however, have been discovered by radial velocity 
measurements alone; none of them were discovered via the 
transit method, nor have any been later discovered to tran- 
sit their host stars. The study of multiple planet systems 
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enables greater understanding of theories of planet forma- 
tion and migration, and affords us the opportunity to study 
planetary dynamics in action, as well as helping us answer 
fundamental questions such as 'how common is the Solar 
System?'. 

Planets that transit their host stars allow us to measure 
fundamental properties such as the planetary radius and re- 
move much of the uncertainty on the value of the planetary 
mass by constraining the orbital inclination angle. This is 
just as true for multiple-planet systems, and further prop- 
erties such as dynamical evo lution time-scal es can be mea- 
sured for transiting systems (|Fabrvckvll2009l ). 
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1.2 Detecting transiting multiple planet systems 

It has been noted that, given that there are now in excess of 
50 known transiting planets, there is a good chance that at 
least one of these systems may har bour additional planets 
which we should be able to detect (|Fabrvckvll2009l ). There 
are three methods for detecting f urther planets in known 
transiting systems (|Fabrvckvll200ij ) . namely (i) searching for 
transit timing variations (TTV) or variations in other tran- 
sit parameters; (ii) searching for long-term radial velocity 
trends; and (iii) searching for additional transits. 

In a multi-planet system, one planet can have a per- 
turbing effect upon the orbit of a second planet, the ef- 
fects of which can include small variations in the timings 
of transits caused by the second planet, s uch that the tran- 
sits are no longer spac ed periodically (|Agol et al ] l2005l : 
iHolman fc MurrEivll2005l ). Many searches for further planets 
have been conducted, and continue to be conducted, using 
TTV, which has great sensitivity to planets in resonant or- 
bits with the first planet, even if the second planet has an 
extremely low mass. 

It is also possible to infer the presence of additional 
planets by measuring lo ng-term trends i n othe r transit pa- 
rameters. For instance, ICoughlin et al.l (|2008t ) report ob- 
served increases in the orbital inclination, transit width and 
transit depth of Gl 436 b, which may indicate the presence 
of another planet. 

In general, the discovery lightcurves of transiting plan- 
ets, such as those produced by SuperWASP, are of insuffi- 
cient quality to measure transit timings and other parame- 
ters with the required precision to discover additional plan- 
ets; predicted timing variations, for instance, are typically 
on the order of seconds or tens-of-seconds. Additional re- 
sources must therefore be expended to obtain high precision 
lightcurves. 

Secondly, known planets continue to be monitored for 
long-term trends in the radial velocity data, as is common 
practice with planets discovered by that means alone. Many 
longer-period additional planets have been found around 
stars around which a relatively short-period planet has 
been discovered by radial velocity measurements. Long-term 
monitoring of RV systems has yielded planets with periods 
of several years (the longest such period is 14.3 yr). This 
too, requires a modest expenditure of telescope time in or- 
der to obtain radial velocity data points at suitable intervals 
to detect longer-period planets. 

Finally, photometric monitoring of known transiting 
systems may reveal transits caused by a second planet in 
the system. The inclination angle of this second planet must 
be sufficiently close to that of the first planet, that the sec- 
ond planet can be seen to transit its star as well as the first. 
This method, unlike the other two, does not necessarily re- 
quire the allocation of further observing resources; instead 
the data archives of transit surveys can be searched for such 
transit signals. Such surveys often observe the same fields 
for several seasons, and so have a large quantity of data on 
known transiting systems. In this paper we present results 
from a search of the data archive of SuperWASP (Wide An- 

5 1e Search for Plane ts), one such wide-field transit survey 
Pollacco et~aLll2006l ). 



1.3 Detecting transiting multiple planet systems 
by transit photometry 

Only two of the 58 known (to 2009 March) transiting planets 
have orbital periods, P, greater than 10 d (these were both 
detected initially by radial velocity means and were only 
later discovered to transit); indeed only six transiting plan- 
ets have periods longer than 5 d. This is largely due to the 
selection effects present in wide-field surveys: (i) in general, a 
relatively large number of transits are required to boost the 
signal-to-noise sufficiently that the trans it may be detecte d 
in the presence of correlated ('red') noise l|Smith et al.ll2006l ). 
This requirement for many transits leads to the preferen- 
tial detection of short-period planets that exhibit frequent 
transits, (ii) The probability that a given planetary system 
exhibits transits is inversely proportional to the orbital semi- 
major axis, a, and so this again causes a bias towards the 
detection of short period planets. 

Any additional transiting planet is likely to orbit with a 
period greater than the currently known planet. Several fac- 
tors militate against the usual difficulties in detecting long- 
period (P > 5 d) planets, however. Most significantly, the 
probability that any second planet will transit should be 
greatly enhanced by the fact that there is already one tran- 
siting planet in the system. This is because the orbits of 
exoplanets in multiple systems are generally predicted to be 
close to coplanar, because such systems are believed to form 
from a fiat disc in a similar fashion to the Solar System (e.g. 
[Bouwman et al. 200(i). In the Solar System, this results in 
all planetary orbits being coplanar to within a few degrees. 

If we assume co-planarity to within 5°, the probabil- 
ity that the second planet transits is approximately equal 
to the ratio of the semi-major axes of the inner and outer 
planets, ain/flout (jPabrvckv 2009!). This means, for example, 
the probability of a planet orbiting a solar analogue with a 
10 day period is increased from about 5 per cent to around 
45 per cent if an inner transiting planet exists with a period 
of 3 days. 

Some previous attempts have been made to detect ad- 
ditional transiting planets. Croll et al. (2007a,b) used the 
MOST (Microvariability and Oscillations of STars) satellite 
to place upper limits on the presence of additional transiting 
planets in two systems (HD 209458 and HD 189733) known 
to harbour a transiting exoplanet. They were able to rule 
out the presence of additional planets larger than about 0.2 
Rj orbiting with periods less than 14 d for HD 209458 and 
planets larger than about 0.15 Pj orbiting with periods less 
than 7 d for HD 189733. 

This paper reports results from an extensive search of 
archival SuperWASP data for additional transits of stars 
known to host transiting exoplanets. We provide constraints 
on the existence of additional transiting planets in such sys- 
tems. 



2 OBSERVATIONS 

Time-series photometry of 24 stars which host transiting 
exoplanets was obtained by the SuperWASP instruments, 
which are wide-field, mu l ti-cam era survey instruments de- 
scribed in iPoUacco et al] (|2006D . 

Fifteen of these planets were observed by SuperWASP- 
N, located at the Roque de los Muchachos Observatory on 
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La Palma in the Canary Islands, and nine by WASP-South, 
at the South African Astronomical Observatory in Suther- 
land. The objects observed are the first eighteen planetary 
systems discovered by SuperWASP (with the exception of 
WASP-9, for which follow-up observations are still ongo- 
ing) and seven similar systems discovered by other transiting 
planet surveys and retrospectively detected in SuperWASP 
data. Details of all these objects are given in Table [1] Su- 
perWASP has an extensive archive of data on these objects. 
Many of them have been monitored for several observing 
seasons, dating back to 2004 in some cases. 



3 SEARCH FOR ADDITIONAL PLANETS 

To search the SuperWASP lightcurves for additional transits 
we first take all existing lightcurves on a particular object 
from the SuperWASP archive. Several lightcurves may exist 
if the star has been observed in several seasons and/or with 
more than one camera. 

A systematic re-analysis of each of these systems was 
performed using all the photometric data either publicly 
available or held by the WASP consortium. The tran- 
sits caused by the known planet were modelled using the 
M arkov-chain Monte Car l o (MC MC) technique described 
by ICollier Cameron et al.l (l2007f ) , which uses the analytic 
model of lMandel fc Agoll (j2002l ') (see Table [l] for details of 
the parameters fitted). The resulting fitted parameters are 
within the joint error bars of the best published parame- 
ters. The modelled transits were then subtracted from the 
lightcurves, and the resulting residual lightcurves concate- 
nated into a single residual lightcurve for each object. 

We then search each of these residual lightcurves using 
huntIstar, a modified version of the adapted Box Least- 
Squares (BLS) algorithm used for routine SuperWASP 
trans it hunting l|Collier Cameron et al.l [20061 : iKovacs et al.l 
l2002h . The huntIstar routine searches the lightcurve of 
a single object for transits, and it is able to handle large 
quantities of photometric data spanning multiple seasons 
and cameras, with a finely-sampled period grid. 

The lightcurves of a total of 24 transiting planet host 
stars (see Table [T] for the details of these objects) were 
searched for additional transits with periods of between 5 
and 25 days using huntIstar. This lower limit of 5 d was 
chosen because, with the exception of WASP-8 b with a 
period of 8.16 d, all of the original planets have periods be- 
tween 0.9 and 5.0 d, and so have been searched at sub-5d 
periods previously. The upper limit of 25 d was chosen after 
consideration of the length of the observing baselines and be- 
cause for this period, the orbital inclination angle, i, must be 
greater than about 88° in order for transits to be exhibited 
by a Jupiter-sized planet (Fig. [l]). The probability of such a 
planet transiting is around 2-3 per cent assuming nothing 
about the system, and about 25 per cent if co-planarity to 
within 5 degrees with a P = 3 d planet is assumed. 

A periodogram is produced for each object, and the 
parameters of the five strongest peaks in the periodogram 
are refined. The resulting candidates are loosely filtered ac- 
cording to the criteria usually used f or SuperWASP planet 
hunting (|Collier Cameron et al.|[2006l ). These criteria are (i) 
At least two transits must be detected; (ii) the reduced 
of the model must be less than 2.5; (iii) the phase-folded 
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Figure 1. Plot indicating the maximum period for a planet to 
exhibit transits, as a function of orbital inclination angle for a 
range of stellar masses. The host star is assumed to be on the 
main sequence and the planet is assumed to be Jupiter-sized. 

lightcurve must not consist of a high proportion of gaps; 
and (iv) the signal-to-red-noise ratio (the S^ed statistic de- 
fined bv ICollier Cameron et al] \200dt )) must be -5 or less. 
We do not apply more strict filtering since we have a small 
number of objects, so the risk of false positives is small. 



4 RESULTS OF SEARCH FOR ADDITIONAL 
PLANETS 

Of the five best peaks for each of the 24 planets (120 peaks in 
total), 75 pass the criteria outlined in Sec. |3l but our initial 
rejection criteria lean strongly towards retaining candidates. 

As we are dealing with a small number of objects, we 
choose not to apply further cuts on, for example signal- 
to-noise, but instead visually inspect all 24 periodograms 
and all 120 phase-folded lightcurves. The periodograms pro- 
duced for each of the 24 stars were inspected for any strong 
peaks indicative of a genuine transit, and the phase-folded 
lightcurves were checked for transit-like signals. Most of the 
objects do not display any noteworthy periodogram peaks; 
Fig. [2] is a typical such periodogram, whereas Fig. |3] is a 
periodogram typical of a genuine transiting planet. None of 
the phase-folded lightcurves display any transit-like signal. 

Five phase-folded lightcurves were also plotted for each 
object, corresponding to each of the five strongest peaks in 
the periodogram. These were inspected for signs of a transit 
at phase 0, but no such signals were observed, even in the 
handful of objects that exhibit at least one reasonably strong 
periodogram peak. 

4.1 WASP-10 

Although a strongish peak at about 12 d is observed in the 
periodogram for WASP-10 (Fig. no transit is seen in 
the corresponding phase-folded lightcurve. This peak is less 
well-defined than the typical peak produced by planetary 
transits (Fig. [5}. Instead, we attribute the peak to stellar 
rotati on, the period, Prot, of which is known to be about 
12 d ^Christian et al.ll2009l ). We confirm the stellar rotation 
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Table 1. Planetary systems searched for additional transiting bodies. The instrument used to observe the system is indicated in column 
2: 'N' represents SuperWASP-N on La Palma and 'S' represents WASP-South in South Africa. Columns 3 to 10 are the parameters of 
the known planet used to subtract the transits of the first planet. They are the orbital period, P, the epoch of mid-transit, to, the transit 
duration, ui, the impact parameter, b, the depth of the transit, (ijp / )^, the stellar density, , the orbital inclination angle, i, and 
the planetary radius, Rp . Also shown is the number of data points in the lightcurve, Wphot- 



Star Obs. MCMC fitted parameters used to subtract transits from lightcurves 
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Figure 2. Periodogram output of HUNTlSTAR for the residual 
lightcurve of WASP-1, after the subtraction of transits due to 
WASP-1 b. This is typical of the periodograms which exhibit no 
strong peaks. 
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Figure 3. Periodogram output of HUNTlSTAR for the unadul- 
terated SuperWASP lightcurve of WASP-1. This is typical of the 
periodogram indicating the presence of a transiting planet, in this 
case WASP-1 b, which orbits with a period of about 2.52 d. 



hypothesis by fitting a sine curve of the form 5+Asm{(jjt + 9), 
where uj = 27r/Prot to the data (Fig. [5]). The best-fitting 
parameters for data taken in SuperWASP-N's 2004 and 2006 
observing seasons are given in Table (2] We find the same 
period of rotation (11.95 d) in both seasons of data, but both 



the phase of rotation and the amplitude of the variability 
differ between the tw o seasons. We also comp ute the auto- 
correlation function (|Edelson fc KrolikI [l988l ) of the data; 
the period determined by this method is 11.84 d. 

Using a generalised Lomb-Scargle periodogram, as de- 
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Figure 4. Periodogram output of huntIstar for WASP-10. The 
strong peak observed at about 12 d is caused by stellar rotation 
with that period. 
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Figure 5. Stellar rotation of WASP-10. The 2004 data are phased 
on a period of 11.949984 d (upper panel) and the 2006 data on a 
period of 11.946772 d (lower panel). Overplotted in each case is 
the best- fitting model of the form <5 -|- As'm{u!t + 9); the values of 
5, A and are given in Table [2] 



Table 2. Best fitting parameters of a sine curve fitted to the 
lightcurve of WASP-10. The parameters are described in the text. 
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scribed in lZechmeister fc Kiirsterl ()2009h . we are able to cal- 
culate the false alarm probability for the signal detected by 
sine fitting. Despite our lightcurve consisting of 8546 points, 
because of red noise the effective number of independent 
data points, rioff, is significantly smaller. To calculate ncff, 
the individual nights of data are shuffled, destroying the ro- 
tation signal, but maintaining the red noise. By assuming 
the highest peak in the resulting periodogram has a prob- 
ability of 0.5, we calculate n^^ = 949. The number of in- 
dependent frequencies i s calculated to b e 1135, according 
to the approximation of Gumming (12004). Using these val- 
ues, and equation 24 of 



Zechmeister fc Kiirsteil (|2009l ). we 



calculate the false alarm probability to be 3.8 x 10"^'^. 

This extremely low probability confirms the reality of 
this variation, as does the fact that very similar periods 
were detected in different seasons of data, and with different 
methods. Taking the mean of all the periods detected, we 
conclude that the star rotates with Prot = 11.91 ±0.05. This 
variability is likely to be caused by starspots; such variability 
is not unusual amongst K-dwarfs. 



5 MONTE CARLO SIMULATIONS 

In conclusion, none of the 24 lightcurves we analysed show 
any evidence for further transiting planets. Although no ad- 
ditional planets were detected in Sec.U it is a useful exercise 
to quantify our ability to detect such planets, and hence to 
determine upper limits to the sizes and orbits of planets we 
have ruled out. To do this, we test our ability to detect plan- 
etary transits with various parameters, through the use of 
Monte Carlo simulations. 



5.1 Generation of artificial lightcurves 

The two parameters that most affect the ability of a survey 
such as SuperWASP to detect a planet are the size of the 
planet and the period of its orbit. We therefore choose to 
determine our ability to detect additional planets as a func- 
tion of these two parameters, following an approach similar 
to that used by the MOST team to place upper limits on 
the presence of additi onal companions i n the HP 209458 a nd 
HD 189733 systems (|CroU et al.ll2007l a: ICroll et al.l l2007^). 

We take the residual lightcurves described in section 
[3] and inject artificial tra nsits into them, using the small 
planet appr oximation o f Mande l fc Agoll l|2002l ). The m- 
puts to the iMandel fc A gol ( 2002) model are the stellar 
mass, stellar radius, effective stellar temperature (Tcff), stel- 
lar limb-darkening parameters and the planet radius and 
orbital inclination angle. We use values of M« , J?* , and 
Tcff taken from the exoplanet.eu website and the non-linear 
limb-darkening coefficients of ^Claret (2000), and we adopt 
an orbital inclination of 89°. Model planetary transits are 
injected at a range of 15 different orbital periods and 10 
different planetary radii, giving a grid of 150 models (see 
Table [3] for all values of P and 7?p ). One hundred different 
lightcurves are created for each of these points in P — J?p 
space, each with a randomly generated epoch of mid-transit, 
to- 



5.2 Searching for injected transits 

Each of the 15,000 lightcurves generated for each object is 
searched for transits by hunt 1 STAR, in exactly the same way 
as were the real data (Sec. |3]). The detection efficiency for 
a planet of a given size and orbital period is determined by 
the fraction of transits recovered for planets of those char- 
acteristics. 

Whilst we inspected visually the periodograms and five 
best lightcurves of each of the 24 planets when searching the 
real data (Sec. |4]), this is clearly an unfeasible proposition 
for 15,000 lightcurves per object. Instead, for lightcurves 
identified as candidates by huntIstar, we require that the 
injected period, or an alias thereof, and the injected epoch 
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Table 3. Model planet parameters used in simulations. Ten different planetary radii and fifteen periods give a total of 150 combinations. 



Periods used / d: 

5.100 6.529 7.957 9.386 10.814 12.243 13.671 15.100 16.529 17.957 19.386 20.814 22.243 23.671 25.100 
Radii used / Rj : 

0.40 0.55 0.70 0.85 1.00 1.15 1.30 1.45 1.60 1.75 



of mid-transit are successfully recovered in at least one of 
the five best periodogram peaks. To do this, we define the 
statistic, 

-Pmcas -fini /h\ 

V=- 5 -, (1) 

-Tinj 

where Pmcas and Pi^j are the orbital periods measured by 
huntIstar, and injected into the lightcurves, respectively. 
We then require rj < 0.001 (corresponding to a detec- 
tion at the injected period), or 0.4995 < rj < 0.5005 or 
—0.0005 < r;* < 0.0005, where rj* = nint{ri)—r] (correspond- 
ing to detections at an alias of the injected period). These 
thresholds were designed to encapsulate clearly-defined pop- 
ulations of objects clustered around ?7 = 0, 77 = 0.5 and 
T) — 1,2,.... Furthermore, we require that Ato < 0.10 d, 
where Ato = \toinj — tomeas\, and toinj and tomeas are the 
injected and measured epochs of mid-transit, respectively. 

Periodograms and phase-folded lightcurves for several of 
the injected transits detected in this manner were inspected 
in the same fashion as the real data (Sec. |3]), in order to 
ensure that these objects could have been detected without 
prior knowledge of the orbital period. The periodogram and 
recovered lightcurve of one such injected transit is shown in 
Fig.H 

5.3 Results of simulations 

5.3.1 WASP-1 

The results of our simulations of extra planets in the WASP- 
1 system are presented in a series of 2 dimensional cuts 
through the parameter space (Figs.[7]&[8ll. Fig.[7]shows de- 
tection efficiency plotted against the radius of the simulated 
additional planets for a variety of orbital periods, whereas in 
Fig. [S] detection efficiency is shown as a function of orbital 
period for a range of orbital periods. 

The full dataset is shown as a contour map in Fig. |5J 
although the results from the models with periods of 7.957 
and 17.957 d are excluded as they are very close to an inte- 
ger number of days. The sharp drops in detection efficiency 
observed at these periods (Fig. [8]) are manifestations of the 
well-known 1 d alias phenomenon, which can cause either 
a sharp reduction or sharp increase in detection ability at 
peri ods which are almo st exactly an integer number of days 
(e.g. ISmith et al.ll2006l ). 

5.3.2 Other systems 

Similar simulations were conducted for several of the 24 sys- 
tems which were searched for additional transits in Sec. |31 
Contour plots of the same type as Fig.[9]were also produced 
for these systems, although not all of them are shown here. 



Periodogram: WASP-4 




-0.1 -0.05 0.05 0.1 

Orbital phase 

Figure 6. An example of an artificial transit which is successfully 
recovered in our simulations. Shown are the HUNTlSTAR results 
from transits corresponding to a planet with P = 15.100 d and 
-Rp = 0.85i?j , which were injected into the WASP-4 lightcurve. 
The planet is detected with P = 15.10138 d in the strongest 
peak of the periodogram (upper panel). The lightcurve, folded 
on the recovered period using the recovered epoch of mid-transit, 
exhibits a clear transit at phase = (lower panel). 



for reasons of space. Our ability to detect additional planets 
in these other systems is very similar to that for WASP- 
1; there are variations in detection efficiency, but this cor- 
relates with the length of the original lightcurve. There is 
more data on WASP-1 than most of the other objects - 
13,630 data points spanning the 2004, 2006, and 2007 ob- 
serving seasons. At the other extreme is HAT-P-4, which 
has a significantly shorter SuperWASP lightcurve, compris- 
ing just 6,546 data points, a few hundred of which are from 
2006, with the rest from 2007. The contour plot for HAT-P-4 
is shown in Fig. 1101 where it can be seen that the detection 
efficiency is poorer, particularly at relatively long periods, 
than for WASP-1 (Fig.[9l). 
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WASP— 1 : Detections vs. radius: various orbital periods 
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Figure 7. WASP-1 simulation results (i). Detection efficiency as 
a function of planetary radius for a second planet orbiting WASP- 
1 for several orbital periods. 
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Figure 10. HAT-P-4 simulation results. Contour map showing 
detection efficiency as a function of orbital period and planetary 
radius. 
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Figure 8. WASP-1 simulation results (ii). Detection efficiency as 
a function of orbital period for a second planet orbiting WASP-1 
for several planetary radii. 
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Figure 9. WASP-1 simulation results (iii). Contour map showing 
detection efficiency as a function of orbital period and planetary 
radius. 



6 DISCUSSION 

6.1 Limits of simulations 

When considering the reliability of the upper limits estab- 
lished in Sec. 15.31 an obvious question to ask is 'could sys- 
tems in which multiple planets exhibit transits have been de- 
tected by SuperWASP (and HAT, TrES, and XO) in the first 
place?'. Although objects which exhibit photometric vari- 
ability are generally rejected as planet candidates, we argue 
that systems with a sub-5d transiting planet and a further 
transiting planet with a period greater than 5 d would be 
unlikely to be rejected. In general, very few transits caused 
by the longer period planet would be present in the discov- 
ery lightcurve and that, combined with the shallow depth 
of planetary transits would prevent the candidate from fail- 
ing test designed to eliminate variable stars. This is borne 
out by simulations conducted on lightcurves containing the 
transits caused by two planets, where the huntIstar algo- 
rithm was still able to detect the inner planet, despite the 
presence of additional transits. 

In particular, we take the example of additional tran- 
sits corresponding to a large, fairly short-period planet {P = 
12.243 d; Rp = 1.45i?j ) which had been successfully recov- 
ered by huntIstar when injected into the residual WASP-4 
lightcurve (Sec. [5]). These same artificially generated tran- 
sits were injected into the unmodified WASP-4 b lightcurve, 
in order to ascertain whether WASP-4 b could still be de- 
tected. WASP-4 b was indeed still detected, with only a very 
slightly reduced signal-to-red-noise value (-13.390 compared 
to -13.564). 

6.2 Lack of detections 

It has been suggested that the lack of transiting multiple 
planet systems is perhaps surprisin g given the prev alence 
of RV multiple systems. Specifically, iFabrvckvl ([200^) notes 
that 11 companion planets are known to exist in the systems 
with the 33 shortest period planets. An apparently similar 
sample of transiting planets, however, co ntains no systems 
with known additional planets. Although IFabrvckvl (|2009l ) 
does acknowledge the existence of observational biases which 
could partly explain this lack of transiting multiple planet 
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systems, we argue that a careful examination of the char- 
acteristics of the known RV multiple planet systems reveals 
that the dearth of such systems is not surprising. 

Of the 28 known multiple planet systems, only four 
have an inner planet which might reasonably have been de- 
tected by a wide-field transit survey, i.e. has P < 5 d and 
Mp sini > 0.2 Mj , suggesting a radius large enough to be de- 
tected. None of these systems contains further planets with 
periods conducive to detection by means of transits; the sec- 
ond planet from the star orbits with P > 95 d in each of the 
four casefl. In other words, there is currently no multiple- 
planet system of the kind we have demonstrated we are able 
to detect at present with SuperWASP known by RV stud- 
ies. This does not, of course, mean that such systems do not 
exist, or that we should not look for them; especially given 
the low expenditure and potentially high reward involved. 



7 CONCLUSIONS 

We conducted a search for additional planets, with periods 
between 5 and 25 d, orbiting 24 stars known to harbour 
a transiting hot Jupiter, using SuperWASP photometry to 
search for additional transits. No planets were detected, so in 
order to place upper limits on the existence of such planets, 
we performed Monte Carlo simulations of planets of various 
sizes, and with various orbital periods. The results of these 
simulations suggest that, for objects like WASP-1 with three 
seasons of SuperWASP photometry, we have a good chance 
(> 50 per cent) of detecting Saturn-sized (Rs = 0.843 Rj ) 
planets out to about 10 d, and a sporting chance of detecting 
such planets with longer periods (there is a ~ 20 per cent 
chance of detecting a Saturn analogue in a 20 d orbit). These 
detection thresholds improve with increasing planet radius 
(up to about 1.2 Rj ), and are lower for stars, like HAT-P-4, 
with fewer data. 

We are able to detect planets larger than about Saturn 
size, and with periods up to ~ 20 d with reasonable effi- 
ciency. As expected, the main factor affecting detectability 
is the time span of the lightcurve used to search for addi- 
tional transits. 

7.1 Future prospects 

As our simulations (Sec. 15. 3p demonstrated, the longer the 
span of the lightcurve, the greater our sensitivity to longer- 
period planets. As transit surveys continue to observe some 
of the known planets, our ability to detect additional tran- 
sits in these systems will increase. Another possibility for 
increasing the number of available photometric data points 
is to share da t a bet ween transit surveys, as suggested by 
[Fleming et all l|2008h . 
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